Diffractive quarkonium production in association with a photon at the LHC 
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The rates for diffractive quarkonium plus prompt-photon associated production at the LHC are 
estimated. The theoretical framework considered in the J/$ (or T) production associated with a 
direct photon at the hadron collider is the non-relativistic QCD (NRQCD) factorization formalism. 
The corresponding single diffractive cross section is computed based on the hard diffractive scattering 
factorization supplemented by absorptive corrections. Such processes are sensitive to the gluon 
content of the Pomeron at small Bjorken-x and they may also be a good place to test the different 
available mechanisms for quarkonium production at hadron colliders. 

PACS numbers: 13.60.Hb, 12.38.Bx, 12.40.Nn, 13.85.Ni, 14.40.Gx 



I. INTRODUCTION 

In recent years, a wide variety of small- a; and hard 
diffractive processes which are under intense study ex- 
perimentally has been used to improve our knowledge 
about Quantum Chromodynamics (QCD). An outstand- 
ing process is the diffractive heavy quarkonium produc- 
tion, where their large masses provide a natural hard 
scale that allows the application of perturbative QCD. 
There are several mechanisms proposed for the quarko- 
nium production in hadron colliders [l], Q , as the color 
singlet model, the color octet model and the color evapo- 
ration model. An important feature of these perturbative 
QCD models is that the cross section for quarkonium 
production is expressed in terms of the product of two 
gluon densities at large energies. In a similar way, the 
diffractive quarkonium production is also sensitive to the 
gluon content of the Pomeron at small- a; and may be par- 
ticularly useful in studying the different mechanisms for 
quarkonium production. From the experimental point of 
view, the heavy quarkonium production is of special sig- 
nificance because they have an extremely clean signature 
through their leptonic decay modes. 

As long diffractive processes are concerned, in hadron 
collisions they are well described, with respect to the 
overall cross-sections, by Regge theory in terms of the 
exchange of a Pomeron with vacuum quantum numbers 
[H]. However, the nature of the Pomeron and its reaction 
mechanisms are not completely known. A good chan- 
nel is the use of hard scattering to resolve the quark 
and gluon content in the Pomeron Such a parton 
structure is natural in a modern QCD approach to the 
strongly interacting Pomeron. The systematic observa- 
tions of diffractive deep inelastic scattering (DDIS) at 
HERA have increased the knowledge about the QCD 
Pomeron, providing us with the diffractive distributions 
of singlet quarks and gluons in Pomeron as well as the 
diffractive structure function In hadronic collisions, 
we shall characterize an event as single diffractive if one 
of the colliding hadrons emits a Pomeron that scatters 
off the other hadron. Hard diffractive events with a large 



momentum transfer are also characterized by the absence 
of hadronic energy in certain angular regions of the final 
state phase space (rapidity gaps). The events fulfilling 
the conditions of large rapidity gaps and a highly excited 
hadron remnant are named single diffractive in contrast 
to those in which both colliding hadrons remain intact as 
they each emit a Pomeron. the so called central diffrac- 
tive events. 

Here, we focus on the following single diffractive pro- 
cesses pp — > p + (J/'fy + 'y) + X and pp — > p+ (T + 7) + X 
at the LHC energies and predict their diffractive ratios 
as a function of transverse momentum (jpr) of quarko- 
nium state. Such processes are interesting because the 
produced large px quarkonia are relatively easy to detect 
through their leptonic decay modes and their transverse 
momenta are balanced by the associated high energy pho- 
ton. As we will verify, the ratio of single diffractive pro- 
duction cross section to the inclusive production cross 
section is not strongly sensitive to the heavy quarkonium 
production mechanism. This channel has been recently 
investigated in the context of heavy ion collisions Q and 
it can also to be compared to the previous calculations 
on diffractive quarkonium production [7j. 

We start by the hard diffractive factorization, where 
the diffractive cross section is the convolution of diffrac- 
tive parton distribution functions and the correspond- 
ing diffractive coefficient functions in similar way as the 
inclusive case. However, at high energies there are im- 
portant contributions from unitarization effects to the 
single-Pomeron exchange cross section. These absorp- 
tive or unitarity corrections cause the suppression of 
any large rapidity gap process, except elastic scattering. 
In the black disk limit the absorptive corrections may 
completely terminate those processes. This partially oc- 
curs in (anti)proton-proton collisions, where unitarity is 
nearly saturated at small impact parameters Q. The 
multi-Pomeron contributions depends, in general, on the 
particular hard process and it is called survival proba- 
bility factor. At the Tevatron energy, s/s = 1.8 TeV, 
the suppression for single diffractive processes is of or- 
der 0.05-0.2 [1, [IS [II El, EDI, whereas for LHC energy, 
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J^s = 14 TeV, the suppression appears to be 0.06-0.1 
[a, E3, E3, El ■ Therefore, these corrections are quite im- 
portant for the reliability of predictions for hard diffrac- 
tive processes. 

This Letter is organized as follows. In next section 
we summarize the main formulas considered to com- 
pute the diffractive ratios for the hadroproduction of 
J/4>(orT) + 7 final state at the LHC energy. In the 
last section we present the numerical results for the in- 
clusive e diffractive cross sections as a function of trans- 
verse momentum at the central rapidity region and give 
predictions for the corresponding diffractive ratios. We 
discuss in details the theoretical uncertainties involved 
in the present estimations and also compare them with 
previous calculations in the literature. 



II. QUARKONIUM PLUS PROMPT-PHOTON 
ASSOCIATED HADROPRODUCTION 



Let us introduce the main formulas for the inclusive 
and single diffractive differential cross sections for the 
production of J/^f + 7 state in proton-proton collisions 
at high energies. For the inclusive case, we will rely on 
the non-relativistic QCD (NRQCD) factorization formal- 
ism and for the diffractive process the calculation are 
based on the Ingelman-Schlein (IS) model for diffrac- 
tive hard scattering fij]. Accordingly, we will take into 
account absorption effects by multiplying the diffrac- 
tive cross section by a gap survival probability factor. 
Concerning inclusive production, as long as J/ip + 7 is 
produced at small longitudinal momentum fraction, the 
gluon fusion channel dominates over the qq annihilation 
process. The corresponding signal is the production of a 
J I ip and an isolated photon produced back-to-back, with 
their transverse momenta balanced. The leading order 
(LO) cross section is obtained by convoluting the par- 
tonic cross section with the parton distribution function 
(PDF), g{x, hf)i in the proton, where fj,p is the factoriza- 
tion scale. At NLO expansion on a s , there are one virtual 
correction and three real corrections processes, as shown 
in Ref. [HIGH- In the NRQCD factorization formalism, 
non-perturbative aspects of quarkonium production are 
organized in an expansion in powers of v, which is the 
relative velocity of the qq in quarkonia. The partonic 
level subprocesses for associated J/ip + j production are 
composed of the gluon fusion subprocesses, which are 



g + g - 7+(cc)[ 3 S 1 \ 3 S 1 8 ] 
9 + 9 - 7 +(cc)[ 1 S 8 , 3 P 8 ] 



(1) 



Quarks initiated subprocesses are strongly suppresses 
and then neglected further. 

In order to obtain the transverse momentum (pt) dis- 
tribution for the process g + g — ► J/ip + 1 + 9, we express 
the differential cross section as 



d 2 cr inc f 2 2 Axix 2 pT da 

— = J dx 1 g p {x 1 ^ F )g p {x 2 ^ F ) 2xi __ Tey Ti 



where we have defined xt = 2mr/-\Aj with yfs be- 
ing the center of mass energy of the pp system and 

nix = \Jpt + rn^ being the transverse mass of outgoing 

J/ip. The gluon distribution, g p (x,Q 2 ), in the proton is 
evaluated at factorization scale fip ■ The common trans- 
verse momentum of the outgoing particles is px and y 
is the rapidity of outgoing J/ip having mass m^. The 
variables x\ and x 2 are the momentum fractions of the 
partons, where M 2 / s < x\ < 1 (M is the invariant mass 
of J/ip + 7 system) and X2 can be written in terms of 
other variables as 



X\Xt& v — 2t 
x 2 = , with r 

2xi — xt£ v 



(3) 



In the NRQCD formalism, the cross section for the 
production of a quarkonium state H is written as <j{H) = 
Yin c n(0|O^|0), where the short-distance coefficients c„ 
are computable in perturbation theory. The (Q|O^|0) are 
matrix elements of NRQCD operators of the form 

(0|O*|0) = £5>|4|iT(A) + X){H(X) + A>„|0). 

X A 

(4) 

The bilinear in heavy quark fields which cre- 

ates a QQ pair in a state with definite color and an- 
gular momentum quantum numbers. Hereafter, we will 
use a shorthand notation in which the matrix elements 
are given as (0% g s ( 2S+1 Lj)). The angular momentum 

quantum numbers of the QQ produced in the short- 
distance process are given in standard spectroscopic no- 
tation, and the subscript refers to the color configuration 
of the QQ: 1 for a color singlet and 8 for a color octet. 
For example, the parton level differential cross sections 
relevant for hadroproduction of Jhp + 7 for the color- 
singlet contributions is given by [l6|, G3] : 
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-s( + t 2 t\ + u 2 u\ 
s\t\u\ 



(Of^ 3 ^)) 



,(5) 



and the full expressions for the color octet contributions 
can be found, for instance, on Refs. [16L Il7|. 

In Eq. ([5]) we have defined 00 = Tr 2 e 2 aa 2 m c /s 2 (with 
charm quark mass m c — 1.5 GeV) and si = S — Am 2 , 
t\ = t — Am 2 , and u\ — u — Am 2 . In these formulas, 
s, t, and u are the Mandelstam variables, which can be 
written as: 



xix 2 s, t = m 



X2\T sm T^ y , u = m,0 — xi^mxe y , 



,(2) 



In our numerical calculations the one loop expression 
for the running coupling, a s (/ii?), with Aqcd = 0.2 
GeV and rif = 4 is considered. The (renormalization 
and factorization) scale for the strong coupling and for 
the evaluation of PDFs is y? F = ii 2 R = {p\ + m|)/4. 
For numerical values of the NRQCD matrix elements 
we have used those from Ref. [l8j, which are (units of 

GeV 3 ): (0^*( 3 5i)) = 1.16, (0^( 3 S X )) = 1.19 x lO" 2 , 



(O^So)} = (O^'^P^/ml = 0.01. We have 
checked that using another set of color octet matrix ele- 
ments, taken from [19], our results do not change con- 
siderably. The production of T(15) -I- 7 can be ob- 
tained from the expressions (J2]{5]) above, by replacing the 
charm mass and charge by the bottom ones (to;, = 4.7 
GeV, eb = —1/3), the J /if) mass by the T(15) mass, 
and by using the corresponding matrix elements. We 
use the values taken from [2(3], namely (units of GeV 3 ): 
(OU 3 Si)) = 10.9, (Oj( 3 Si)) = 0.02, {OjfSo)) = 0.136 
and (OJ( 1 P )) = 0. 

Lets move now to the hard diffractive cross section. 
The associated J/ip + 7 single diffractive production at 
large pr consists of three steps: first a hard Pomeron is 
emitted from one of the protons in a small squared four- 
momentum transfer \t\. After that, partons from the 
Pomeron interact with partons from the other hadron. 
Finally, J/tp + 7 (or T + 7) are produced in the final 
state, from the point-like QQ by the soft gluon radia- 
tion. Considering the px distribution and based on the 
IS model for diffractive hard scattering, the differential 
cross section can be expressed as 
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FIG. 1: The inclusive (upper solid line) and single diffractive 
(bottom lines) cross section for J/SP + 7 hadroproduction as 
a function of pr at the LHC. The diffractive cross section is 
computed for two sets of diffractive gluon distribution (FIT 
A and FIT B). Only central rapidity region, \y\ < 1, is con- 
sidered. 
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dx\ I dt /:p/p(:Enp, t) 



2 2 Ax 1 x jp x 2 Pt da 

x gsp{x l ,^ F )g p {x 2 ,H F )- = — r-Tf,6 

2x\Xsp — xre y dt 



where Xtp is the momentum fraction of the proton car- 
ried by the Pomeron, t is the squared of the proton's 
four-momentum transfer, x\ is the momentum fraction of 
the partons inner the Pomeron and x 2 is the momentum 
of the partons inner the anti-proton. Here, fw/p(xw,t) 
is the Pomeron flux factor, which can be written as 

fw/p( x TP,t) oc x]p 2a ^F 2 (t). In the following calcula- 
tion, we use the Pomeron factor flux considered by the 
HI Collaboration Q analysis of the diffractive structure 
function in diffractive DIS process. The variable x 2 and 
the corresponding Mandelstam variable read now as 
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2x\Xtp — XT& y 
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(7) 

(8) 
(9) 



To calculate the hard diffractive cross sections, Eq. 
((6j> , we consider a standard Pomeron flux from Regge 
phenomenology and which is constrained from the ex- 
perimental analysis of the diffractive structure function 
|5J], where xtp dependence is parametrized as 



2a P (t)-l 



(10) 



with the Pomeron trajectory assumed to be linear, 
ap(t) = ajp(0) + a' w t and Bjp, a' w and their uncer- 
tainties are obtained from fits to HI FPS data 12111. 



A-p is the normalization parameter, choosed such that 
xtp ■ J" t ™* n fw/pdt = 1 at Xtp = 0.003, where \t min \ w 
TOpXp/(l — xjp) is the minimum kinematically accessible 
value of \t\, m v is the proton mass and \t cut \ = 1.0 GeV 2 



is the proton mass and \t cut \ 
is the limit of the measurement. For the diffractive gluon 
distribution in the Pomeron, gsp{x\, ^ 2 F ), we will consider 
the diffractive PDFs obtained by the HI Collaboration 
at DESY-HERA @, where the Pomeron structure func- 
tion has been modeled in terms of a light flavor singlet 
distribution £(#), i. e., the u, d and s quarks with their 
respective antiquarks. Also, it has a gluon distribution 
g(z), with z being the longitudinal momentum fraction of 
the parton entering the hard subprocess with respect to 
the diffractive exchange. The gluon density is a simple 
constant at the starting scale for evolution, which was 
chosen to be Q\ = 2.5 GeV 2 . In our numerical calcula- 
tions, we will use the following cuts for the integration 



< Xtp < 0.05, where x 



over Xtp, ip"" < Xtp < 0.05, where i™ n — 

As a final ingredient for the single diffractive cross 
section, we also consider the suppression of the hard 
diffractive cross section by multiple-Pomeron scattering 
effects (absorptive corrections). This is taken into ac- 
count through a gap survival probability < \S\ 2 >, which 
can be described in terms of screening or absorptive cor- 
rections [22I ]. Let A(s, b) be the amplitude of the partic- 
ular diffractive process of interest and considering it in 
the impact parameter space b, the probability that there 
is no extra inelastic interaction is 



< |5| 2 >= 



fd 2 b\A(s,b)\ 2 exp[~n(s,b)} 
fd 2 b\A(s,b)\ 2 



(11) 



where fi is the opacity (or optical density) of the inter- 
action of the incoming hadrons. There are intense the- 
oretical inves tig ation on this subject in last years. We 
quote Ref. [131 ] for a nice discussion and comparison 



of theoretical estimations for the gap survival probabil- 
ities. We notice that it is the main theoretical uncer- 
tainty in the present calculation of diffractive ratios. As 
a baseline value, we follow Ref. [12], which considers 
a two-channel eikonal model that embodies pion-loop in- 
sertions in the Pomeron trajectory and high mass diffrac- 
tive dissociation. For LHC energy on pp collisions, one 
has < \S\ 2 >= 0.06. In single channel eikonal models, 
this factor can reach up to 0.081 — 0.086 as discussed in 
Ref. [Hj]. Concerning the model dependence, the single 
channel eikonal model considers only elastic rescatter- 
ings, whereas the multi channel one takes into account 
also inelastic diffractive intermediate re-scatterings. The 
available experimental observables which can be com- 
pared to the theoretical predictions of the survival prob- 
ability factor are the hard LRG di-jets data obtained in 
the Tevatron and HERA [l2|, HH as well as diffractive 
hadroproduction of heavy bosons (W and Z°) in the 
Tevatron [13|. 



A. Results and comments 
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FIG. 2: The inclusive (upper solid line) and single diffractive 
(bottom lines) cross section for T (15) +7 hadroproduction 
as a function of pr at the LHC. The diffractive cross sec- 
tion is computed for two sets of diffractive gluon distribution 
(FIT A and FIT B). Only central rapidity region, |y| < 1, is 
considered. 



Let us now present the estimates for the quarkonium 
plus prompt-photon associated diffractive hadroproduc- 
tion for the LHC energy. In Fig. [T]we present the numeri- 
cal calculations for the inclusive cross section (upper solid 
line) for J/ty + 7 production as a function of transverse 
momentum, given by Eq. (|2|), using the MRST gluon 
PDF [24J. The numerical results include both the color 
singlet and color octet contributions to the production 
process. The corresponding branching ratio into dilcp- 
tons is taken into account and we concentrate on the cen- 
tral rapidity region, imposing the kinematic cut \y\ < 1. 
The order of magnitude is not small, giving units of nano- 
bars in the fully integrated case. Notice that the absolute 
value of inclusive cross section is strongly dependent on 
the quark mass, NRQCD matrix elements and on the 
factorization scale and it is know pjl HH that the NLO 
corrections strongly enhances the cross section at large 
Pt- This theoretical uncertainty is minimized when we 
compute the corresponding diffractive ratio. The single 
diffractive cross section (without absorption corrections) 
is shown in the bottom lines. We present the estimate 
using two different sets of the diffractive gluon distribu- 
tion function determined experimentally. The result for 
FIT A parameterization is represented by the dot-dashed 
curve, whereas the FIT B calculation is given by the thin 
dashed curve. It is verified that the absolute value of 
diffractive cross section is weakly sensitive to the uncer- 
tainties on the diffractive gluon PDF, being more pro- 
nounced in the charmonium case and on the small pt re- 
gion. We perform a similar analysis for production of the 
bottomonium state T(15). In order to do so, in Fig. [5] 
the inclusive cross section and the single diffractive cross 
section for production of final state are presented as a 
function of pr- The label for the curves follows those al- 
ready considered in Fig. [1] The cross sections are about 



one order of magnitude smaller than the charmonium 
case. This is understood from the dependences on quark 
charge and mass a cx Cq/toq, which gives roughly speak- 
ing a^/ay cx (el / el)[ml / m 2 ] = 4[m^/m^]. As already 
mentioned, the single diffractive cross section is almost 
insensitive to the diffractive gluon PDF in the range on 
Pt presented here. The situation would be different for 
the central diffractive case, i.e. the double Pomeron ex- 
change (DPE) process, which is dependent on the prod- 
uct of diffractive gluon distribution on both Pomerons. 
Notice that our numerical results are somewhat consis- 
tent with those presented in Ref. [25[ , where the J/'I' + 7 
cross section is computed in the same framework and a 
renormalized Pomeron flux is considered for the investi- 
gation of single diffractive cross section. Our estimate 
for the Upsilon (diffractive) case has been not addressed 
before in literature. 

Our main results are presented in Table I, where the 
prediction for the diffractive ratio as a function of px 
is shown (using FIT A for the diffractive PDFs). As a 
sample of numerical results for the inclusive and single 
diffractive cross section, we also give the estimates for 
charmonium and bottomonium plus prompt-photon as- 
sociated hadroproduction in the LHC regime. The trans- 
verse momenta are presented in the range 4 < px < 10 
GeV for central rapidity region \y\ < 1. It was verified 
that when considering the forward rapidity region (for 
instance, \y\ > 2) a smaller diffractive ratio is found. 
This is consistent with previous calculation, presented 
in Ref. (25|. However, we have found a slightly large 
diffractive ratio in comparison to that work. The reason 
can be twofold: for the inclusive case, we are using a dif- 
ferent factorization scale (fip = Et in Ref. (2f|) and in 
the diffractive case we considered the eikonal-inspired ab- 
sorptive corrections against renormalized Pomeron flux 
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TABLE I: Inclusive and single diffractive cross sections in 
units of picobarns for the interval 4 < px < 10 GeV for central 
rapidity region (integrated over \y\ < 1). The corresponding 
diffractive ratios, Rsd, are presented considering the KKMR 
value for the absorptive corrections at the LHC energy, < 
|Sf >= 0.06. 



p T [GeV] 
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10 




97.04 


36.46 


14.54 


6.21 


2.82 


1.36 


0.68 


^ (m 


0.78 


0.26 


0.10 


0.04 


0.017 


0.008 


0.0036 


Rsd [%} 


0.8 


0.71 


0.69 


0.64 


0.6 


0.59 


0.53 




5.91 


3.88 


2.49 


1.58 


1.00 


0.64 


0.41 


^ (T) 


0.036 


0.022 


0.013 


0.008 


0.0054 


0.003 


0.0018 


#SD [%] 


0.6 


0.56 


0.53 


0.51 


0.54 


0.47 


0.44 



As a summary, the theoretical predictions for inclu- 
sive and single diffractive J/ip + 7 and T + 7 produc- 
tion has been done for the LHC energy in pp collisions. 
The estimates for the differential cross sections as a func- 
tion of quarkonium transverse momentum are presented, 
discussing the main theoretical uncertainties. The corre- 
sponding diffractive ratio is computed using hard diffrac- 
tive factorization and absorptive corrections. It is verified 
that the ratios are less dependent on the heavy quarko- 
nium production mechanism and quite sensitive to the 
absolute value of absorptive corrections. We found that 
at the LHC, the diffractive ratio for J/tp+j state is about 
Rsd = 0.8 - 0.5 % in the interval 4 < p T < 10 GeV and 
-Rsd = 0.6 — 0.4% in same interval for the T + 7 state. 



[26j taken in Ref. [25| (the authors also neglected the Q 2 
evolution in the gluon density of the Pomeron). This last 
fact could also to explain the pt dependence appearing 
in our estimates for the diffractive ratios, whereas in Ref. 
[25| such a dependence is weaker. 
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